Simple analysis of off-axis solenoid fields using the scalar magnetostatic potential: 
application to a Zeeman-slower for cold atoms 
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In a region free of currents, magnetostatics can be described by the Laplace equation of a scalar 
magnetic potential, and one can apply the same methods commonly used in electrostatics. Here we 
show how to calculate the general vector field inside a real (finite) solenoid, using only the magnitude 
of the field along the symmetry axis. Our method does not require integration or knowledge of the 
current distribution, and is presented through practical examples, including a non-uniform finite 
solenoid used to produce cold atomic beams via laser cooling. These examples allow educators to 
discuss the non-trivial calculation of fields off-axis using concepts familiar to most students, while 
offering the opportunity to introduce important advancements of current modern research. 

PACS numbers: 03.75.Be, 03.75.Lm, 84.40.Az, 73.21. Cd 



I. INTRODUCTION 

Magnetic fields produced by solenoids and axially sym- 
metric coils are ubiquitous, and the ability to calculate 
them is an integral part of training in physics. Time 
constraints, however, tend to focus the attention of most 
introductory electromagnetism (EM) courses to the an- 
alytical solution of only a few highly symmetrical cases, 
such as the field along the axis of a circular coil or inside 
an infinite solenoicP"^]. Nevertheless, many applications 
require at least an estimate of the full vector field in re- 
gions away from the axis^', which involve mathematical 
tools often not discussed at the introductory level. On 
the other hand, most EM courses already dedicate a fair 
amount of time teaching students to identify and solve 
electrostatic problems using the Laplace equation. In 
some cases the same methods can be applied to magne- 
tostatic problems, sometimes leading to useful insights. 

Sadly, most students do not usually appreciate the 
similarities between the two classes of problems^ due 
to a limited exposure to practical examples involving 
the magnetostatic potential. We feel that this ability 
is usefuP, particularly because scalar potentials are gen- 
erally more intuitive and easier to visualize. Besides, a 
unified treatment could be pedagogically relevant ingen- 
cralizing the discussion of the multipole expansions"-^. 
Therefore, the primary goal here is to present a couple of 
pedagogical examples illustrating the application of the 
magnetostatic potential method to real solenoids. 

In addition, these examples also offer the opportunity 
to discuss in the classroom axisymmetric fields evaluated 
off-axis, without the need to introduce the formalism of 
elliptic integrals. Although other methods for finding 
off-axis magnetic fields have been mentioned earlier in 



the literature^ 13 , to our knowledge this has not been pre- 
sented from such a simple and intuitive viewpoint. 

Moreover, as further motivation, we have chosen an ex- 
ample that brings a real and practical application from 
the cutting edge of research into the classroom: a non- 
uniform solenoid used in many research laboratories to 
produce beams of slow (cold) atoms. This solenoid, called 
Zeeman-sloweJ^K^, is used in conjunction with appro- 
priately prepared laser beams to slow down and cool neu- 
tral atoms, from hundreds of Kelvin to milliKelvin tem- 
peratures, by combining the action of radiation pressure 
with the Zeeman effect. This device is one o f the staple 
developments in the area of laser coolingS^ 3 , and one of 
the enabling technologies leading to the 1997 Nobel prize 
in PhysicsIS The techniques for laser cooling and trap- 
ping of atoms have produced many dramatic advance- 
ments in our understanding of quantum physics^, in- 
cluding the achievement of Bose-Einstein condensation^, 
which was recognized with another Nobel prized in 2001. 
In both cases, magnetic fields have been an important 
part of experimental design and data interpretation. Ed- 
ucators can astutely use the solenoid discussed here, as 
well as the references herein, to introduce and discuss 
some of these modern developments in quantum physics, 
making the subject even more interesting to students. 



II. REVIEWING SOME BASIC CONCEPTS 

We begin here by recalling the fundamental equation 
of magnetostatics : , where H is the mag- 

netic field and J the current density. Typically if is 
related to the magnetic induction field if by some con- 



stitutive relation expressing the properties of a particu- 
lar material. For linear and isotropic materials, with a 
magnetic permeability /i, B = and in a current-free 
region = 0, implying that B = —v(f>M- Since 

Maxwell's equations also state that V • B = 0, this re- 
sults in S7 2 cf>M = 0, which is Laplace's equation for the 
magnetic potential (j>M, m any current-free region. 

Although Laplace's equation is only typically valid in 
a region free of charges or currents, these are allowed to 
exist on or outside a surface S surrounding that region. 
The solutions of Laplace's equation present three impor- 
tant properties: superposition, smoothness and unique- 
ness. The property of superposition results from the fact 
that Laplace's equation is a linear equation. Smoothness 
implies that no solution in a region V of space, bounded 
by a surface S, can present either a maximum or a min- 
imum within V (extreme values can occur only at the 
surface S). The thirdproperty is the one most relevant 
to us here, as it states^ that if one finds a solution 4>m, in 
a region of space consistent with the prescribed bound- 
ary conditions, that solution is unique up to an additive 
constant. Therefore, it does not matter what particular 
method is used to find the solution. Once an appropriate 
solution is found, it is uniquely valid. 

However, despite the obvious similarities between the 
electrostatic and magnetostatic potentials, there ar e in- 
deed reasons why the analogy can only be taken so faiM^I 
and is not widely explored further in textbooks. The first 
one arises whenever J ^ 0, in which case it is not trivial 
to write a relation between <pM and J . The second com- 
plication occurs due to the fact that the scalar potential 
is generally a multiply valued function, requiring a pre- 
scription specifying where it can be used. However, as it 
has been shown by Bronzan^, these complications can 
be overcome, permitting one to exploit the advantages of 
the concept of a scalar magnetic potential. 



III. THE MAGNETIC FIELD OF A FINITE 
UNIFORM SOLENOID 

We start by considering the field in the interior of a fi- 
nite uniform solenoid carrying a current /, as illustrated 
in Fig. [T] For generality and convenience, we describe 
the problem using spherical coordinates. In this geome- 
try it is easy to note that due to the axial symmetry of 
the problem the field B(z) depends only on z — rcos9. 
The magnetostatic potential can be found using for a 
boundary condition the magnitude of the field along z, 
which is readily available through simple summation for- 
mulas over the approximately circular coils forming the 
solenoid. 

In spherical coordinates, the solution of the axisym- 
metric scalar potential 4>m can be written in the form: 
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Figure 1: Schematic representation of a finite and homoge- 
nous solenoid. The relevant coordinates are presented in (a), 
while their details are shown in (b). The crosses (dots) rep- 
resent current flowing into (out of) the plane of the page. 

where a£ and bi are coefficients to be determined and the 
Pi represents a Legendre polynomial of order I. 

Because we are mainly interested in the values of the 
field inside the solenoid, we set bg — to avoid the sin- 
gularity at r = 0. As a result, the potential takes the 
simpler form: 

4>M{r,9) = a + a 1 rP 1 (cos 9) + a 2 r 2 P 2 (cos 9) + .... (2) 

For points along the z-axis, we have cos 9 = 1 and 
Eq. ([2| becomes 

4>m{z) = a + aiz + a 2 z 2 + a 3 z 3 + (3) 

Incidentally, we can in general also expand the scalar 
potential in a Taylor series around some point Zq, 



<I>m{z) = 4>m(zq)+(z-z )- 



dz 



(z - z ) 2 d 2 <j) M 
2! 



dz 2 |Z0 " 

and comparing it to Eq. ([3]), for zq = 0, we obtain the 
coefficients ag in terms of the series expansion 



aa = 0m(O); 



1 d e <j> M (z) 
i\ dz 1 



z=0 ■ 



(5) 



In this way, the full scalar potential becomes analytically 
determinable, allowing us to evaluate ~B* = —v4>m a t 
any point in space where Eq. ^ applies. 

As a first example, let us now consider the case of 
a finite solenoid of length L and radius R, carrying a 
uniform current /. If the solenoid has N turns per unit 
length, the magnetic field along the z-axis can be easily 
calculated by integrating the expression for the axial field 
of a circular current loop^, resulting in: 



3 



where z± = z ± |r, and a = /iqNI/4tt in SI units. Now, 



B z (z) = - 



>M 



dz 



we can write 



<j> M {z) = - / B(z')dz'. 



(7) 



(8) 



Using Eq. ^ in Eq. ([8]) we obtain the general form of 
the potential for the finite solenoid along the axis: 



4>m(z) 



[R 2 + z\) 2 - (R 2 + z 2 



(9) 



Expanding Eq. ^ around z — 0, as in Eq. @, we 
get the various coefficients for (j) M (z). Using these coef- 
ficients and introducing the expression for the Legendre 
polynomials Pi(cos9), while keeping terms up to third 
order, we finally get: 
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From the last equation, one can calculate the compo- 
nents Bg and B r of the magnetic field by simply taking 
the gradient of the potential: 
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dr 
aL cos 9 



(11) 
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It is interesting to note that these results give approxi- 
mate analytical results for the magnetic field everywhere 
inside the solenoid, with their precision limited by the 
number of terms included in the power series expansion. 
One can test these results by comparing the expressions 
(111 and ( 12 ) with those presented in Chapter 5 of refP, 
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Figure 2: Numerical calculation of the exact solution (using 
elliptic integrals) for the field of a finite solenoid of length 
L = 1 m, radius R — 10 cm, number of windings N — 100 
and carrying a current I = 1 A. The points (black) show the 
magnitude |B| along the axis (p = 0), whereas the solid (blue) 
line represents B z , and dashed (red) line B p at p = 8 cm. 



where a different method was used to evaluate the field 
components. In particular, we will show that if one keeps 
only the first order in the expansion, the result simplifies 
to the approximate solution of problem 5.2 of reference^. 
For that we recall the relations 



Bp = B r sin 



Bg cos t 
B z = B r cos 9 — Ba sin t 



(13) 



from which we obtain, up to third order, 



B„ 



3aLR 2 r 2 
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(14) 



Finally, using rain 9 — p and rcos9 = z, in the limit 
R/L < 1, Eq. (Tub yields 



B p (z,p) 



96ttNIR 2 pz 
c U' 



(15) 



which is expressed here in CGS units, with a — 2nNI/c, 
to facilitate direct comparison with the result presented 
in refP. 

To show a practical application of the method, we now 
compare our results to a realistic numerical calculation of 
a finite solenoid. Figure [2] shows the numerical results for 
the magnitude of the magnetic field and its components. 
The calculation assumes that the solenoid is composed 
of a series of single circular coils, and performs a direct 
summation over the exact analytical expression, based on 
elliptic integrals, for each individual coil. The numerical 
result was determined to be very close to the experimen- 
tal results, as at this scale the error introduced by ap- 
proximating the actual helical winding by a sequence of 
circular single coils is typically negligible. The relevant 
physical parameters are given in the captions of Fig(2] 
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Figure 3: Comparison of numerical (points) and analytical 
(lines) results for the uniform solenoid. In (a) the axial and 
in (b) the transverse field profiles have been shown. All the 
results were evaluated off-axis, at p = 8 cm, and the analytical 
results are shown for expansions up to third (dashed-black) 
and fifteenth (solid-red) orders. The inset in (a) shows the 
good agreement obtained near z=0. 



Further, in Fig. |3j we compare the analytical results 
obtained by keeping the first eight terms in the expan- 
sion (corresponding to the fifteenth order in r) against 
the numerical calculations shown in Fig. [2j In addition, 
Fig. [3] also shows the good partial agreement obtained 
using the third order approximation, extending to dis- 
tances up to about half the size of the solenoid. Note 
that, for a real finite system, the disagreement increases 
rapidly after some point. That can be improved signif- 
icantly by including higher order terms, allowing for a 
much better approximation near the edges, as shown in 
Fig. [3] However, due to some of the simplifications made, 
the approximate analytical result still does not contain 
all the physics of the problem. For instance, it does not 
accurately describe the field outside the solenoid. Never- 
theless, the magnetostatic potential method still provides 
a reasonable representation of the internal fields up to the 
very end of the solenoid. 
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Figure 4: (a) Sketch of a tapered (triangular shape) solenoid 
creating an inhomogeneous current distribution to produce 
the appropriate axial field profile (b) for a Zeeman-slower. 



IV. THE ZEEMAN-SLOWER: AN 
INHOMOGENEOUS FINITE SOLENOID 

Now we will consider another very interesting and 
practical problem, familiar to many atomic physics labo- 
ratories, namely the design of a solenoid capable of pro- 
ducing an axial field with a parabolic profile, as in 



B(z) = B b + B Q y/l -Pz, 



(16) 



where B b , Bq and j3 are constants. Such a field is suitable 
for slowing atomic beams using laser lighlP. The field 
of Eq. ( 16 ) causes a spatially varying Zeeman effect that 



compensates for the changing Doppler shift of the moving 
atoms, thus keeping them in resonance with the light as 
they decelerate along the beam path. This technique is 
called Zeeman slowing^, and the parabolic shape is cho- 
sen to keep the radiation pressure constant, typically at 
a rate of ~ 10 6 raj s 2 , throughout the Zeeman solenoicP^ 
shown in Figure |4j 

In general, the atomic beam encompasses a certain 
solid angle as it traverses the solenoid and most atoms 
follow trajectories which do not lie exactly on the axis. 
Since the resonance condition with the laser depends on 
both the magnitude (via the detuning) and direction (via 
the polarization) of the magnetic field, the knowledge 
of the off-axis field is important in understanding how 
light interacts with atoms at different points inside the 
solenoid. 

The magnetic potential along the z-axis, in this case, 
takes the form 

Mz) = ~J B(z')dz' = -B b z+^(l-(3z)i, (17) 
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Figure 5: Numerical results representing the experimental 
field of a real Zeeman solenoid, with L = 125 cm, R = 4 cm, 
and B max ~ 120 mT. The solid (blue) line shows \B\ along 
the axis (p = 0), whereas the dashed (red) line represents the 
transverse field B p at p = 2.5 cm. 



where the constant of integration has been suppressed. 



Following the same steps in section III and after cal- 
culating the derivatives and solving for the coefficients 
a^, we obtain the general form of the magnetic potential 
for the Zeeman solenoid: 
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where T(n) is the gamma function. Now we can calculate 
the spherical components of the magnetic field, 



B r (r,6) = {B b + B ) cose 
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x [cos6 l P„(cos6') — P n _i(cos9)] . 
The transverse and axial components can be easily ob- 



tained from Eq. (13), with the shorthand z — 
B p (p, z) 
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Notice that the transverse component of the magnetic 
field inside the Zeeman-slower does not depend on B^. 
Also, note that B„ = at p = (on-axis), as expected. 
It can be verified that the on-axis field sums back to the 
exact expression of Eq. (16). Although caution may be 
necessary when evaluating the field for p = (9 = 0), the 
careful use of L'Hospital's rule ensures finite and correct 
answers. 

Now, to compare these analytical approximations to 
the experimental field shown in Fig. [5j we will follow a 
different approach here. The motivation is to mimic a 
situation where the current distribution that generates 
the field may not be known exactly, but the axial field 
can be measured directly in the laboratory. This could 
be the case in a real practical application, where imper- 
fections in the winding pattern often are not considered 
in the ideal model. From the experimental data one can 
then build a mathematical model, using cither a fitting 
function (if the functional form is known or can be easily 
guessed), or by using an interpolating function, such as 
a polynomial, to represent the data in a limited region 
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Figure 6: Comparison between numerical (from Figure [5| 
and analytical solutions for the Zeeman solenoid for (a) B z for 
p — 0, and (b) B p for p = 2.5 cm. Points (blue)represent the 
experimental field and solid lines (red, purple) the analytical 
approximations of Eq.(22) and Eq.(21) respectively. 
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of space. Here, since the approximate functional form of 
the axial field is known, we will extract the model param- 
eters by numerically fitting the data in Fig.[5j to B z (z) in 
Eq. (16 1, and substituting them in the expressions (21) 
and (22). Note that the limitations of this type of mod- 
eling may result in some inaccuracies, particularly close 
to the edges, where fringe effects are important. In any 
practical situation one may need to explore different ap- 
proaches to find a mathematical model accurate enough 
in the region of interest. 

After following these steps to model the data, we show 
in Fig.|6]a comparison between our analytical approxima- 
tions and the experimental field. There is a reasonable 
agreement between the solid lines, representing equations 
(21) and (22) and the data points, representing the ex- 
perimentally determined field. Note that, in contrast to 
the uniform finite solenoid where the power series was 
used to approximate the exact solution, here the power 
series simply approximates our modepil function. There- 
fore, increasing the order^l of the expansion only im- 
proves the agreement with the model (fitting) function, 
which represents the data only over a limited region and 
still does not contain all the information about the fields 
in the problem. This is clearly visible in Fig. [6] where 
good agreement is found only in the range z ss (0.4 to 
1.2) m. 



CONCLUSION 



Using the simple concept of the magnetostatic scalar 
potential, and only the knowledge of the field along the 
symmetry axis, we have shown how to determine the 
vector magnetic field anywhere inside an inhomogeneous 
finite solenoid, without explicitly integrating (or even 
knowing) the current distribution. In cases where the 
current distribution is known, but the expression for the 
field off-axis is non-trivial (for instance, given by ellipti- 
cal integrals), one can still gain some insight by using the 
method described here. This simple analysis follows from 
a straightforward analogy with the electrostatic bound- 
ary value problem, and can be useful in determining field 
inhomogeneities in various practical experiments involv- 
ing solenoids. In the present article we have used an 
example from contemporary atomic physics experiments 
to demonstrate the method. However, we believe that 
a simplified version of this discussion (e.g.: the uniform 
finite solenoid) could be used in a regular classroom set- 
ting, as a practical example of a calculation of off-axis 
magnetic fields, for undergraduate courses and teaching 
laboratories in physics and engineering. 
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